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Abstract: We investigate the effect induced by neutrino oscillation on the dark matter indirect 
detection signal which consists in a muon neutrino flux produced by neutralino annihilation in the 
Earth core. We consider the neutrino oscillation parameters relevant to the atmosferic neutrino deficit, 
both in the v., — > and v., — > i^, cases. 



in 
cn 

o 



X 



1. Introduction 

Many different techniques have been proposed 
for the detection of dark matter particles which 
could make up the halo of our Galaxy. Among 
the different possibilities, the detection of a neu- 
trino flux by means of neutrino telescopes rep- 
resents certainly an interesting tool, which is al- 
ready at the level of imposing some (mild) con- 
straint on the particle physics properties of the 
neutralino, the most interesting and studied dark 
matter candidate This particle is present in 
all the supersymmetric extensions of the stan- 
dard model as a linear combination of the super- 
partners of the neutral gauge and higgs fields. In 
the present paper we will perform our calcula- 
tions in the minimal supersymmetric extension 
of the standard model (MSSM), for a definition 
of which we refer to Ref . 1^ and to the references 
therein quoted. 



2. Up— going muons from neutralino 
annihilation in the Earth 

The neutrino flux has origin from neutralino pair 
annihilation inside the Earth where these dark 
matter particles can be accumulated after having 
been captured by gravitational trapping. The 



differential flux is 
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where Ta denotes the annihilation rate, d is the 
distance of the detector from the source (i.e. the 
center of the Earth), F lists the neutralino pair 
annihilation final states, B^^ denotes the branch- 
ing ratios into heavy quarks, r lepton and gluons 
in the channel F. The spectra dN ^^^/dEi, are 

/I'm 

the differential distributions of the (anti) neutri- 
nos generated by the r and by hadronization of 
quarks and gluons and the subsequent semilep- 
tonic decays of the produced hadrons. For de- 
tails, see for instance Rcfs. Here we 
only recall that the annihilation rate depends, 
through its relation with the capture rate of neu- 
tralinos in the Earth, on some astrophysical pa- 
rameters, the most relevant of which is the local 
density pi. 

The neutrino flux is produced in the inner 
part of the Earth |Q and propagates toward a 
detector where it can be detected as a flux of 
up-going muons, as a consequence of neutrino- 
muon conversion inside the rock that surrounds 
the detector. A double differential muon flux can 
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Figure 1: Muon response function dN^/d\ogx vs. 
tiie parent neutrino fractional energy x — Ei^/m^ 
for neutralino annihilation in the Earth. Different 
curves refer to different neutralino masses : = 50 
GeV (solid), = 80 GeV (dotted), = 120 GeV 
(shot-dashed), — 200 GeV (long-dash), — 
500 GeV (dot-dashed). 
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Figure 2: Flux of up-going muons <J>J^=''^"^ from neu- 
tralino annihilation in the Earth, plotted as a func- 
tion of m^. The solid line denotes the present up- 
per limit 1^]. Different neutralino compositions are 
shown with different symbols: crosses for gauginos, 
open circles for higgsinos and dots for mixed neu- 
tralinos. 



be defined as 
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where Na is the Avogadro's number, g{E^, E'^; X) 
is the survival probability that a muon of initial 
energy E'^ will have a final energy i?^ after prop- 
agating along a distance X inside the rock and 

^ daAEy,E') 
S{E., E,) = ^ $.(K) (2.3) 

where i = v^,v^ and da^_-){E^, E'^) /dE'^ is the 

charged current cross-section for the production 
of a muon of energy E'^ from a neutrino (antineu- 
trino) of energy E^. 

A useful quantity for the discussion in the 
following Sections is the muon response function 
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where E^^ is minimal energy for detection of up- 
going muons. For SuperKamiokande and MACRO, 
E^^ ~ 1.5 GeV The muon response func- 
tion indicates the neutrino energy range that is 



mostly responsile for the up-going muon signal. 
Fig. 1 shows a few examples of it, plotted as func- 
tions of the variable x = Eu/my^, where my. de- 
notes the neutralino mass. Fig. 1 shows that the 
maximum of the muon reponse happens for neu- 
trino energies of about Ei, ~ (0.4 — 0.6) m^, with 
a half width which extends from E^ ~ 0.1 to^ to 
Ey ~ 0.8 m^. 

Finally, the total flux of up-going muons is 
defined as 



(2.5) 



The natural background for these kind of 
searches is represented by the flux of up-going 
muons originated by the atmospheric neutrino 
flux. Experimentally one searches, inside a small 
angular cone around the center of the Earth, 
for a statistically significant up-going muon ex- 
cess over the muons of atmospheric origin. 
No excess has been found so far and therefore, 
an upper limit on $^ can be derived. Fig. 2 
shows the present most stringent upper limit ob- 
tained by the MACRO CoUaboration [||. In the 
same figure the theoretical calculations of for a 
scan of the supersymmetric parameter space are 
also displayed. The plot refers to pi = 0.3 GeV 
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cm ^ and is obtained by a variation of the MSSM 
parameters in the following ranges: 20 GeV < 
M2 < 500 GeV, 20 GeV < |^| < 500 GeV, 
80 GeV < ruA < 1000 GeV, 100 GeV < mo < 
1000 GeV, -3 < A < +3, 1 < tan/? < 50. 
For further details of the calculation, we refer 
to Ref. [Q. The comparison of the scatter plot 
with the experimental upper limit would imply 
that a fraction of the supersymmetric configura- 
tion could be excluded. However, a variation of 
the value of pi inside its range of uncertainty can 
lower the theoretical prediction by about a factor 
of 3 As a consequence, we can conservatively 
consider that only a small fraction of the susy 
configurations can be potentially in conflict with 
the experimental upper limit, when no oscillation 
effect on the neutrino signal is assumed. 

3. Neutrino oscillation effect on the 
up— going muon signal 

The recent data on the atmospheric neutrino deficit 
indicate that the may oscillate, either into 
Vt or into a sterire neutrino I's |, |7|. If this 
is the case, also the i'^ produced by neutralino 
annihilations would undergo an oscillation pro- 
cess. The energies involved in both atmospheric 
and neutralino-produced neutrinos are the same. 
The baseline of oscillation of the two neutrino 
components is different, since atmospheric neu- 
trinos cross the entire Earth, while neutrinos pro- 
duced by neutralino annihilation travel from the 
central part of the Earth to the detector (we 
recall once more that neutralinos annihilate in 
the core of the Earth). On the basis of the fea- 
tures of the Vfj, oscillation which are required to 
fit the experimental data on atmospheric neutri- 
nos d, 0, we expect that also the neutrino flux 
from dark matter annihilation would be affected. 
In the next Sections we will explicitely discuss 
the i/f^ — > and the —^ Vs cases, in a two 
neutrino mixing scenario 1^. 

3.1 Vr vacuum oscillation 

In the case of — > oscillation, the z/^ flux 
is reduced because of oscillation, but we have 
to take into account also that neutralino anni- 
hilation can produce which in turn can oscil- 
late into i^^ and contribute to the up-going muon 
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Figure 3: i/p survival probability in the case of i/^ — > 
Vr oscillation. The solid line refers to sin^(2S„) = 1, 
the dashed line is for sin^(2S„) = 0.8. In both cases, 
Am^ = 5 ■ 10^=^ eV 

N. Fornengo (1999) 
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Figure 4: Scatter plot of the ratio 'y'^^SttSwI^i^ vs. 
the neutrahno mass m^. ('l'p)J|!3^ni is the up-going 
muon flux in the case of — > oscillation, while 
is the corresponding flux in the case of no oscillation. 

flux. The Vr flux can be calculated as discussed 
in Sect. 1 for the flux, and it turns out to be 
always a relatively small fraction of the flux. 
The muon neutrino flux at the conversion region 
can therefore be expressed as 
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where the vacuum survival probabiUty is 



N. Fornengo (1999) 



l-sin^(26'^)sin^ 



(3.2) 
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where Am^ is the mass square difference of the 
two neutrino mass eigenstates, 9^ is the mixing 
angle in vacuum and R is the Earth's radius. Fig. 
3 shows the survival probability for two differ- 
ent values of the neutrino oscillation parameters. 
Smaller (larger) values of Ato^ have the effect 
of shifting the curves to the left (right). Com- 
paring Fig. 1 with Fig. 3, we notice that the 
reduction of the up-going muon flux is stronger 
when there is matching between the the energy 
El ~ 5.2 • 10"3^m2(cV^) of the first (from the 
right) minimum of the survival probability and 
the energy Ei, ~ 0.5m^ which is responsible for 
most of the muon response in the detector. This 
implies that a maximum reduction of the signal 
could occur for neutralino masses of the order of 
mx(GeV) ~ 104Am2(eV^). The i^r ^ os- 
cillation makes the reduction of the muon flux 
less severe, but it is not able to completely bal- 
ance the reduction effect because the original Vr 
flux at the source is sizeably smaller than the Ur 
flux. Therefore, the overall effect of the neutrino 
oscillation is to reduce the up-going muon sig- 
nal. This effect is summarized in Fig. 4, where 
the ratio between the up-going muon signals in 
the presence and in the absence of oscillation 
are plotted as a function of the neutralino mass. 
The susy parameter space has been varied in the 
same ranges quoted for Fig. 2. We notice that 
the strongest effect is present for light neutrali- 
nos, since in this case the muon flux is mostly 
produced from neutrinos whose energy is in the 
range of maximal suppression for the oscillation 
phenomenon. The effect is between 0.5 and 0.8 
for < 100 GeV. On the contrary, the fluxes 
for larger masses are less affected, and the reduc- 
tion is less than about 20% for > 200 GeV. 



3.2 j/„ 



matter oscillation 



In the case of i/^ oscillation, the neutrino 

flux is simply 
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Figure 5: survival probability in the case of —> 
Vs oscillation, for s\v?{29v) = 0.8 and Am^ — 5-10~'' 
eV ~^ . The solid line refers to neutrinos, the dashed 
line is for antineutrinos. 
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Figure 6: Scatter plot of the ratio (4'n)^rin/'^'M vs. 
the neutralino mass m^. ('l'^i)^^n is the up-going 
muon flux in the case of — > Vs oscillation, while "1>^ 
is the corresponding flux in the case of no oscillation. 



and no i/^ regeneration is possible from the ster- 
ile neutrino. In this case, the effective potential 
of and Vg inside the Earth are different and 
we have to solve the evolution equation for prop- 
agation in the core and in the mantle. Neutrinos 
(produced in the center of the Earth) cross once 
half of the core and once the mantle. By consid- 
ering both core and mantle as of constant density, 
we can express the survival probability as ||, [l0|] 



(3.4) 
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where U is the 2x2 neutrino mixing matrix, 9a 
{a = c,m for core and mantle, respectively) are 
the effective mixing angles in matter and they 
are related to the vacuum mixing angle 9^ as 



sin2(20,) = 



s\n\29,)ea 



[(e,cos(20,) + l)2+e2sin2(20„)] 

(3.5) 

with = AmV(2i;,K); Va = ±Gi.7V;j/V2 
is the matter potential in a medium of num- 
ber density iV," for neutrinos (+) and antincutri- 
nos (— ). In Eq.(3.4), D is the evolution matrix 
Dij{4>a) = <5ydj, where dl = 1, d^ = exp(z0a) 
and 

(ba = VaRa [(^a cos(2(?„) + 1)^ + sm\29,)] 

(3.6) 

In Fig. 5 an example of the z/^ and i^^ sur- 
vival probability is given for representative val- 
ues in the range allowed by the fits on the atmo- 
spheric neutrino data Am^ = 5 • 10~^ eV^ 
and sin^(26'i,) = 0.8. From Fig. 5 and the previ- 
ous discussion relative to Fig. 3, we expect that 



in the case of Va 



the reduction of the muon 



signal is significantly less severe than in the case 
of I'p, — !■ Vt- In fact, in these case the minima 
of the survival probability occur for lower neu- 
trino energies, and threfore the oscillation can 
affect only muon fluxes originated by very light 
neutralinos. This is manifest in Fig. 6, where 
the ratio of the up-going muon fluxes in pres- 
ence and absence of oscillation are shown. In 



oscillation the theoretical predictions are re- 
duced in the presence of oscillation. With the 
oscillation parameters deduced from the fits on 
the atmospheric neutrino data, the effect is al- 
ways larger for lighter neutralinos. In the case 



of Va 



the reduction is between 0.5 and 0.8 



for < 100 GeV and less than about 20% for 
> 200 GeV. In the case of 



this casa, the leducliou of the signal ib alwa^is 



Vs, the re- 
duction of the signal is up to 30% for neutralino 
masses lower that about 80 GeV and smaller than 
10% for heavier neutralinos. 
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